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The average diameter of the first micrometer particle size standard (Standard Reference Material 1690), an 
aqueous suspension of monosized polystyrene spheres with a nominal 1 fj.m diameter, was accurately determined 
by three independent techniques. In one technique the intensity of light scattered by a diluted suspension of 
polystyrene spheres was measured as a function of scattering angle, using a He-Ne laser polarized in the vertical 
direction. The second technique consisted of measuring as a function of angle the intensity of light scattered from 
individual polystyrene spheres suspended in air, using a He-Cd laser with light polarized parallel and perpendic- 
ular to the scattering plane. The measurement of row length by optical microscopy for polystyrene spheres 
arranged in close-packed, two-dimensional hexagonal arrays was the basis of the third technique. The mea- 
surement errors for each technique were quantitatively assessed. For the light scattering experiments, this 
required simulation with numerical experiments. The average diameter determined by each technique agreed 
within 0.5% with the most accurate value being 0.895±0.007 fxm based on light scattering by an aqueous 
suspension. Transmission electron microscopy, flow through electrical sensing zone counter measurements, and 
optical microscopy were also used to obtain more detailed information on the size distribution including the 
standard deviation (0.0095 /xm), fraction of off-size particles, and the fraction of agglomerated doublets (1.5%). 

Key words: array sizing; index of refraction; light scattering; Mie scattering; particle counting; particle size 
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1. Introduction 

This paper describes the measurements performed to 
accurately characterize the average particle size of the 
nominal one-micrometer polystyrene spheres particle 
size standard (Standard Reference Material 1690). The 
spheres were purchased from Dow Chemical Com- 
pany. 1 They have a nearly monosize distribution and are 
dispersed in water at a weight concentration of 0.5%. 

About the Authors: G. W. Mulholland is a research 
chemist in NBS' Center for Fire Research. A. W. 
Hartman, G. G. Hembree, Egon Marx, and T. R. 
Lettieri are physicists in the Bureau's Center for 
Manufacturing Engineering. 

'Certain materials are identified in this paper in order to adequately 
specify the experimental procedure. Such identification does not im- 
ply recommendation or endorsement by the National Bureau of Stan- 
dards, nor does it imply that the materials or equipment identified are 
necessarily the best available for the purpose. 



The three techniques used for measuring average par- 
ticle size were: light scattering from individual spheres 
suspended in air, light scattering from an aqueous sus- 
pension of the spheres, and optical row-length mea- 
surements of spheres arranged in two dimensional ar- 
rays. The selection of these techniques resulted in part 
from a detailed review of all the generic measurement 
techniques appropriate for the dimensional calibration 
of microscopic-particle size standards by Swyt [l] 2 . 
These same three techniques were used previously by 
others in characterizing a batch of monosize polystyrene 
spheres produced by Dow Chemical Co. (batch number 
LS-1028-E). Phillips et al. [2] obtained a diameter of 1.20 
fxm by single particle scattering measurements, Rowell 
et al. [3] obtained 1.21 jam by measurement of light 
scattering from a suspension of particles, and Bierhuizen 
and Ferron [4] obtained 1.20 jam based on array sizing. 
The good agreement among these independent 



2 Figures in brackets refer to literature references at the end of this 
paper. 



measurements gave us confidence that an accurate par- 
ticle size standard was achievable. 

Our study represents an extension of the earlier work 
in regard to the accuracy achievable by the three tech- 
niques. The major effort of our study was directed at 
performing the measurements in such a way as to max- 
imize the accuracy in the determination of the mean 
particle size. This required the use of accurate cali- 
bration standards for magnification and angular align- 
ment in addition to a quantitative error analysis. 

Perhaps the most widely accepted technique for mea- 
suring the particle size of polystyrene spheres has been 
transmission electron microscopy (TEM). This is the 
basis of the particle size information that is provided 
with the monosize polystyrene spheres marketed by the 
manufacturer. We also had intended to use the TEM as 
one of our primary techniques, but we were unable to 
find an accurate method for calibrating the mag- 
nification of the microscope. The relatively wide range 
of values for particle diameter (from 1.10 to 1.26 jam) 
reported by five independent TEM measurements [3] 
for Dow Latex LS-1028-E also suggests that the method 
is unsatisfactory for highly accurate particle size mea- 
surements. While we did not use the TEM for deter- 
mining average size, we did use it for determining the 
standard deviation, cr P , of the size distribution. 

Methods that have been used by others for measuring 
the particle size of polystyrene spheres include flow 
ultramicroscopy, ultracentrifugation, quasi-elastic scat- 
tering, and small-angle x-ray scattering. Some general 
design limitations of these instruments are given below. 
Our choice of light scattering and optical array sizing 
techniques was based on their being at least competitive 
with other methods in regard to sizing accuracy, and on 
our familiarity with the measurement methods and the 
theoretical basis. 

Both the quasi-elastic scattering technique and small- 
angle x-ray scattering work best for particles smaller 
than about 1 jam. The former technique [5] measures, in 
effect, the particle diffusion coefficient, which increases 
with decreasing particle size. The resolution is limited 
for 1 /xm particles but improves with decreasing size. 
This is just the opposite of angular-intensity light scat- 
tering and array sizing where the resolution drops off 
markedly for particle sizes less than about 0.5 jam di- 
ameter. Small-angle x-ray scattering [6] is also best for 
small particle sizes because of the small wavelengths of 
the x-rays. 

Davidson et al. [7] obtained an uncertainty of about 
±4% at the 95% confidence level for flow ultra- 
microscopy based on counting 300 particles. The major 
source of uncertainty was attributed to the statistical 
uncertainty in the average counting rate. The uncer- 
tainty in the number of particles counted is simply equal 



to the square root of that number. Van den Hul and 
Vanderhoff [8] obtained a mean diameter for Dow La- 
tex LS-1028-E about 18% less than the mean of the 
other techniques listed by Rowell et al. [3] using ultra- 
centrifugation. A major drawback in the technique for 
this application is the small density difference, 1.00 ver- 
sus 1.05, between water and the polymer. Nevertheless, 
Van den Hul et al. point out that all the deviation from 
other measurements cannot be accounted for by the 
uncertainty in the density difference. Apparently no de- 
tailed error analysis has been made for the ultra- 
centrifuge applied to particle measurements. 

While the major focus of our study was the accurate 
determination of the number average particle diameter 
D n , we also obtained information regarding the size dis- 
tribution including the standard deviation, o>, of the 
particle size distribution, the fraction of off-size par- 
ticles, and the fraction of agglomerated doublets. In ad- 
dition, we obtained an estimate of the deviation from 
sphericity based on TEM measurements. A brief de- 
scription of the emulsion polymerization technique for 
producing polystyrene spheres and the factors affecting 
the stability of the suspension is given in the Appendix. 



2. Light Scattering From Individual Spheres 

We determined the diameter of the polystyrene 
spheres by measuring the angular distribution of the 
intensity of the light scattered by a single sphere and 
comparing measured intensities with those computed 
for different values of the diameter, D, and the refrac- 
tive index, n . 

The particles were levitated electrostatically in a Dif- 
ferential II light-scattering photometer developed by 
Wyatt and Phillips [9]. This type of instrument has been 
used to determine D and n for spheres with diameters 
1.099 urn (Phillips et al. [2]), 0.796 urn (Cook and Ker- 
ker [10]), 1.011, 0.794, and 0.600 /am (McRae [11]), and 
1.01, 5.7, and 11.9 jam (Davis and Ravindran [12]). The 
sizes given above are nominal sizes provided by Dow 
Chemical Company, the manufacturer of the particles. 
Bottiger et al. [13] have developed the capability of 
measuring all 16 components of the Mueller matrix, the 
light-scattering phase matrix, for a single polystyrene 
sphere. They also collected the particles on the tip of a 
needle for subsequent observation in an electron micro- 
scope. 

A detailed error analysis enabled us to make a quan- 
titative statement regarding the accuracy of this tech- 
nique. This and other aspects of the experiment and data 
analysis have been previously reported by Marx and 
Mulholland [14]. 



2.1 Experimental Method 
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An aerosol made up of polystyrene spheres is gener- 
ated by nebulizing a suspension of the spheres in dis- 
tilled, filtered water; the water evaporates rapidly leav- 
ing charged spheres. One of these spheres is then 
levitated electrostatically in an optical cell with elec- 
trodes as shown in figure 1. The particle is kept in the 
light beam of a He-Cd laser (X=441.6 nm) by means of 
a servocontrol. 

The scattered light reaches a photomultiplier (having 
an S-21 response) via a traveling periscope, which has 
an acceptance angle of about 2° and which moves in a 
horizontal arc of almost 180°. We found that the output 
of the angle encoder deviates by as much as 3° from the 
true angle, and we performed an angle calibration using 
an accurately indexed protractor. 

We use a half-wave retardation plate to change the 
polarization of the incident light. 

A typical sequence of measurements begins with a 
scan of the background, then a particle is captured and 
one or more scans is performed to measure the light 
scattered by that particle. A rotation of the half-wave 
plate by 45° then changes the polarization of the incident 
beam by 90°, and one or more scans is performed to 
measure the scattered light, followed by a background 
measurement for the new polarization. The light 
intensity-angle data pair closest to each integer angle 
over the range 20° to 160° is stored on magnetic tape. 

2.2 Data Analysis 

The data are matched to the results of the calculations 
of the angular distribution of the intensity of the light 
scattered by a uniform sphere, as found by Mie and 
others [15-171. 

The values of D and n are determined by minimizing 
the quality of fit 



Ground 



Figure 1-Single particle scattering 
instrument (Differential II). The 
pillbox shaped scattering cell has 
a pin electrode (Vj) separated by 
an insulator from the ground 
plate electrode and opposite to 
the base plate electrode (V 2 ). 



where N is the number of data points, E t is the measured 
intensity at the angle 6 if T t is the computed value for the 
same angle, and a is an overall scaling factor required 
because we do not measure absolute intensities. The 
expression for the a that minimizes Q is given by 
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The least-squares fit is then determined by a grid 
search in the Dn -plane for progressively finer grids. 
This procedure can be speeded up using an iterative 
search such as the Marquardt-Levenberg algorithm af- 
ter an initial coarse grid search. 

Other measures of the quality of fit are possible, and 
we consider mainly those that involve a change of vari- 
ables. Specifically, we use the transformed intensity 
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because it tends to equalize the height of all the peaks for 
a particle of a size around 1 ju,m in diameter. In the 
original variable 1(6), the large peak at small angles is 
overemphasized. On the other hand, the use of li(0) 
overemphasizes the values at angles close to 180°, which 
are small with large relative errors. 

2.3 Results 

The theoretical best fit curves are plotted in figure 2 
for a single polystyrene sphere with and without the 
weighting function. There are four independent deter- 
minations of the best fit diameter; two based on vertical 
polarization and two based on horizontal polarization 
data. The best fit values range from 0.898 to 0.905 \im, 
which is a range of about 0.8%. A small portion of the 
table of values of Q(n,D) near the minimum is given in 
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Figure 2-Best fits to the intensity of scattered light versus angle data for a single polystyrene sphere as determined by the least squares method 
with (bottom) and without (top) the sfn 2 0/2 weighting factor. Data ( x x); theory ( ). 



table 1. The quantity Q varies rapidly with respect to D 
for fixed n (column) and rapidly with respect to n for 
fixed D (row) but is rather insensitive in the diagonal 
direction along which the product nD is almost con- 
stant. In a three-dimensional plot of \/Q as a function of 
n and D t there is a very sharp peak when the surface is 
viewed in one direction but almost a ridgelike appear- 



ance in the other as shown in figure 3. The relative 
insensitivity of Q to correlated changes in n and D neces- 
sitates that careful angle calibrations and detailed data 
analysis be performed in order to obtain accurate values 
for n and D. 

The results for the eight SRM 1690 particles analyzed 
are contained in table 2 including the averages for the 



Table 1. Normalized values of Q near the minimum for a SRM 1690 particle. 



Vertical Polarization 



Horizontal Polarization 



n 
D,]i.ra 
0.899 
0.901 
0.903 
0.905 
0.907 
0.909 



1.600 

4.75 
2.34 
2.31 

1.65 
L32 
1.32 a 



1.604 

3.18 
2.15 
1,47 
1.10 
1.06 
1.33 



1.608 

2.18 

1.48 

1.09 

I 

1.20 

1.70 



1.612 

1.70 
1.29 
1.15 

1.29 
1.70 
2.40 



1.616 

1.69 
1.50 
1.58 
1.92 

2.52 
3.39 



n 


1.608 


1.612 


1.616 


1.620 


1.62< 


D,[x.m. 












0.893 


14.22 


9.78 


6.36 


4.07 


3.03 


0.895 


9.58 


5.93 


3.37 


2.01 


1.92 


0.897 


5.90 


3.11 


1.51 


1.16 


2.13 


0.899 


3.35 


1.55 


1 


1.81 


4.06 


0.901 


2.15 


1.46 


2.16 


4.35 


8.13 


0.903 


2.59 


3.24 


5.40 


9.25 


14.92 



"The minimum value of Q in each column is underlined to illustrate the slow variation of Q on a diagonal. 




1/Q 



VERTICAL 




HORIZONTAL 



1/Q 




1/Q 



Ridge axis 



Ridge axis 



VERTICAL 




HORIZONTAL 



Figure 3-Two views of the three-dimensional surface of the inverse of the quality of fit, 1/Q, as a function of the index of refraction, n , and the 
diameter, D, over a range of a few percent from the best value. The ridges cause difficulties for an accurate determination of the correct values 
of n and£>. 



Table 2. Radius and index of refraction for eight particles from SRM 1690. 





Harmonic Mean" 




Radius 


Refractive 




urn 


Index 




0.4545 


1.614 




0.4570 


1.613 




0.4490 


1.618 




0.4465 


1.612 




0.4480 


1.613 




0.4505 


1.613 




0.4535 


1.608 




0.4394 


1.606 


Average 


0.4498 


1.612 


cr 


0.0051 


0.004 



Vertical 




Polarization 


Radius 


Refractive 


/im 


Index 


0.4550 


1.614 


0.4570 


1.613 


0.4530 


1.605 


0.4510 


1.600 


0.4515 


1.610 


0.4535 


1.610 


0.4535 


1.608 


0.4402 


1.603 


0.4518 


1.608 


0.0048 


0.005 



Horizontal 


Polarization 


Radius 


Refractive 


/Am 


Index 


0.4535 


1.610 


0.4530 


1.623 


0.4490 


1.618 


0.4460 


1.614 


0.4480 


1.613 


0.4505 


1.613 


0.4535 


1.610 


0.4394 


1.610 


0.4491 


1.614 


0.0045 


0.004 



a Harmonic mean refers to the determination of the best fit radius and refractive index by taking the harmonic mean, Q h , of the Q's for the cases of vertically and 
horizontally polarized light, Q v an d Qh respectively (i.e., 2/Q h =l/Q v +l/Q H ). 



diameter and index of refraction, 0.900 /xm and 1.612, 
respectively. The average particle size in this report 
refers to the number average defined as follows: 



A= 



N 

N ' 



(4) 



where TV refers to the number of particles measured. The 
particle size is not perfectly uniform (o> =0.010) and we 
expect to see a spread in the measured diameters. On the 
other hand, the refractive index should be the same for 
all particles, and we attribute the variation to experi- 
mental error. 



For one particle we took nine scans with vertically 
polarized light and three more with horizontal polar- 
ization. The results of the analysis are shown in table 3. 
The standard deviation for the diameter, 0.0022 jam, is 
smaller than the one previously calculated for groups of 
particles, but that for the index of refraction, 0.005, is 
about the same, suggesting that this parameter does not 
change, at least for particles of the same size. 

2.4 Error Analysis 

To study the effects of experimental errors on the best 
fit parameters, we modified the program that produces 
the theoretical values of the scattered intensities to add 



Table 3. Results for repeated scans for one particle. 



Run 


Unweighted Fit 


Weighted 


Fit 
















n=1.608 a 


n=1.615 




Radius 


Refractive 


Radius 


Refractive 


Radius 


Radius 




jxm 


Index 


/xm 


Index 


/im 


fim 


1-V 


0.4402 


1.601 


0.4380 


1.612 


0.4384 


0.4366 


2-V 


0.4364 


1.620 


0.4376 


1.614 


0.4398 


0.4378 


3-V 


0.4392 


1.606 


0.4366 


1.618 


0.4386 


0.4372 


4-V 


0.4402 


1.603 


0.4398 


1.605 


0.4388 


0.4372 


5-V 


0.4394 


1.606 


0.4410 


1.601 


0.4390 


0.4372 


6-V 


0.4380 


1.610 


0.4398 


1.604 


0.4386 


0.4368 


7-V 


0.4398 


1.604 


0.4386 


1.610 


0.4388 


0.4370 


8-V 


0.4394 


1.606 


0.4388 


1.609 


0.4390 


0.4372 


9-V 


0.4394 


1.607 


0.4400 


1.604 


0.4392 


0.4374 


1-H 


0.4394 


1.610 


0.4402 


1.601 


0.4398 


0.4384 


2-H 


0.4380 


1.612 


0.4402 


1.612 


0.4386 


0.4374 


3-H 


0.4384 


1.614 


0.4364 


1.620 


0.4394 


0.4382 


Average 


0.4390 


1.608 


0.4389 


1.609 


0.4390 


0.4374 


cr 


0.00110 


0.0053 


0.00150 


0.0063 


0.00047 


0.00053 



"The last two columns correspond to values of the radii obtained for a fixed value of the refractive index (unweighted fit). The first value is obtained for the best 
fit to all particles of this size, and the second one is the published bulk value. 



simulated errors and generate a file of "data" to be pro- 
cessed by the same procedure as described above. Noise 
obtained using a random number generator with a max- 
imum amplitude of 1-3% of the intensity at each angle 
was added to all curves. 

The major component of the random error is attrib- 
uted to the noise in the intensity of the scattered light. 
This noise can be monitored by fixing the angle and 
recording a trace of the intensity versus time. The gen- 
eral characteristics of such a trace show that there are 
two components of this noise. The smaller component, 
about 2% of the intensity, affects each angle reading 
independently; the larger component, which we think 
comes from the motion of the particle within the laser 
beam when the servocontrol is activated, is about 5% 
and has a lower frequency so that it affects intensity 
readings in groups of about 10 degrees for the data col- 
lection speed that we used. The uncertainty in particle 
size for a 0.900 jam particle diameter was found to be at 
most 0.3% or 0.0027 jam. As described above, the ob- 
served variability in particle diameter for 12 repeat scans 
of the same particle resulted in a cr of 0.0022 jam, which 
is consistent with most of the error being attributed to 
noise in the intensity of the scattered light. 

The quantity of primary interest is the random error 
associated with the determination of the average par- 
ticle diameter. The standard deviation of the size distri- 
bution, cr P , is much larger than the cr associated with 
repeat measurements and will ultimately limit the accu- 
racy in the mean size determination. Our measured oof 
the particle diameters based on eight particles is 0.010 
jam; an estimate based on our electron microscopy for 
about 160 particles, as discussed below, is slightly 
smaller, 0.0095 jam. The component of the overall un- 
certainty in particle diameter resulting from the width of 
the size distribution is given by 



R =t m (0.025) 






(5) 



where n is the number of particles (eight in this case) and 
t m (0.025) is the Student /-value for m degrees of freedom 
and for 95% confidence level. The number of degrees of 
freedom relates to the determination of a> and is one less 
than the number of particles sized by electron micros- 
copy, 159. We obtain a value of R equal to 0.0066 /xm. 
The sources of systematic errors include angle drift, 
8 g , polarization misalignment, S P , wavelength uncer- 
tainty, d\, and non-volatile impurities, Si. There is an 
error on the order of 2% if the angle encoder output is 
not calibrated. To this is added the problem of a slight 
drift in the calibration that was apparently electronic in 
origin. Over the two days during which the light scat- 
tering measurements were performed on the eight poly- 
styrene spheres, the angle drift was about 0.08°. Includ- 



ing this drift in a numerical simulation, we find the 
uncertainty in particle size due to the drift, 8 e , to be 
0.0028 jam. 

The mixture of polarization states comes primarily 
from imperfections in the optical components and the 
misalignment of the periscope. The imperfections in the 
half-wave plate and other optical components allow the 
leakage of a small fraction of light (1 part in 200) at the 
orthogonal polarization direction. The error from the 
inclination of the periscope with respect to the scatter- 
ing plane is comparable to the error due to the optical 
imperfections. A combined estimate of these effects 
comes to S P =0.0018 jam. 

The least important error is the uncertainty in the 
wavelength of the He-Cd laser, S x = 0.00007 /xm. 

In addition to instrumental errors, there are uncer- 
tainties associated with the nature of the particle. While 
the polystyrene spheres are often considered to be ho- 
mogeneous, we expect the existence of a surface coating 
from the non- volatile emulsifier added during the prepa- 
ration of the particles. If we assume that all the emulsi- 
fier resides on the particle surface, we find that the 
diameter is increased by 0.0008 jam, so that 
di= —0.0008. There is also a contribution to the surface 
coating from the 2 ppm impurities in the distilled water 
in which the particles are suspended, but it is negligible 
in comparison with the emulsifier contribution. 

An estimate of the total uncertainty, U T , is obtained 
by adding the random error, R , to the sum of the abso- 
lute values of the systematic errors. 

tfT=-R + |8*| + |8p| + |8,| (6) 

Using the values given above, we obtain t/ T =0.012 jam 
with a number average size of 0.900 jam. 

The above analysis includes effects for which we 
have quantitative estimates. Two other possible system- 
atic error sources are slight asphericity of the particle 
and possible inhomogeneity of the particle resulting 
from strain as the water evaporates from the particle. 
Transmission electron microscope measurements of the 
particle described in section 5 indicate an asphericity of 
0.6% at a precision of about 0.3%. Over the course of 
our measurements the particle would have undergone 
extensive rotation as a result of Brownian motion. We 
have no quantitative estimate of this effect, though we 
intuitively expect our value to be close to the volume 
equivalent sphere diameter. 

3. Light Scattering From a 
Suspension of Spheres 

Light scattering measurements on a suspension of par- 
ticles is complementary to single-particle scattering 
measurements. The signal-to-noise performance for sin- 



gle-particle measurements is limited by the small scatter- 
ing intensity and by particle movement in the beam. For 
scattering from a particle suspension the signal-to-noise 
ratio is much improved because of the large number of 
particles in the scattering volume and because of the use 
of photon counting. An accurately indexed rotary table 
was used for our light-scattering measurements with a 
suspension, whereas a rapidly-traversing periscope with 
some drift in the angle calibration was used for the 
single-particle measurements. One disadvantage of mea- 
surements with a suspension is that there is one more 
unknown in the measurement, namely the size distribu- 
tion, compared to the single-particle measurements 
where particle size and the index of refraction are the 
only unknowns. 

There have been a number of studies of light scatter- 
ing by nearly monosized suspensions. Kerker's book [17] 
includes several studies of the determination of the mean 
size and standard deviation based on the angular de- 
pendence of the polarization ratio. Our instrument de- 
sign is similar to the one described by Wims and Myers 
[18]. The recent determination of particle size of poly- 
styrene spheres by Rowell et al. [3] by light scattering 
measurements is similar to our study in terms of mea- 
surement technique and data analysis. 

3.1 Design of the Photometer 

A schematic diagram of the light scattering photome- 
ter for particle suspension is shown in figure 4. A 5-mW 
intensity-stabilized He-Ne laser beam is directed to the 
center of a 10-cm diameter scattering cell with optical 



flat windows on opposite sides. Incident light is verti- 
cally polarized and only the vertically polarized com- 
ponent of the scattered light is detected (VV scattering). 
A beam stop is positioned in the cell to minimize the 
amount of reflected light reaching the detector. The 
collection optics (figs. 4 and 5) consist of two apertures, 
a 1 X microscope objective, a vertical polarizer, and a 
glass diffusing screen. The scattered light is detected by 
a low-noise photomultiplier tube with a small (2.5 mm 
diameter) photo cathode. The output pulses from the 
anode are amplified, discriminated, and then counted 
over a 10-second time interval. The typical number of 
counts per 10-second interval at the first peak is about 
200,000, with repeat measurements agreeing within 
about Wo. The lowest count was on the order of 10,000. 

For some measurements, part of the laser beam is split 
off and sent into a reference photomultiplier. Output 
counts from the reference amplifier-discriminator are 
used as clock pulses to the photon counter. In this man- 
ner, a ratio is taken between the signal photon counts 
and reference photon counts, thus cancelling out the 
laser intensity fluctuations. This is not a significant effect 
for the intensity stabilized He-Ne laser, but it is signifi- 
cant for the He-Cd laser to be discussed later. 

A rotary table with an accuracy and reproducibility 
of ±1 arc second (0.0003°) is used for the angle mea- 
surements. The rotary table is indexed at two-degree 
increments. The PMT is mounted rigidly to an arm at- 
tached to the rotary table. Using a micrometer adjust- 
ment and a dial indicator, the center of the scattering 
cell is made coincident with the center of rotation of the 



Figure 4-Light scattering photom- 
eter for particle suspension. Op- 
tical windows not shown. 
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Figure 5-Design of detector optics for minimizing acceptance angle. 
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rotary table. Next the laser is mounted so that the beam 
goes through the center of the two optical flats and 
reflects back on itself. The pinholes in the detector sys- 
tem are then positioned in the center of the laser beam. 
The estimated deviation of the zero angle alignment, 
including a slight error in positioning the cell (about 0.05 
mm) and in the pinhole positioning (about 0.2 mm), leads 
to an uncertainty in the zero angle of about 0.07°. The 
effect of this uncertainty on the particle size deter- 
mination will be discussed in subsection 3.4. 

It is important to have a small acceptance angle for 
the detector in order to obtain good size resolution. 
This is accomplished in our design by using two pin- 
holes before the detector as indicated in figure 5. In 
addition, a collection lens is used to limit the depth of 
field. With these optics in the collection arm, the mea- 
sured acceptance angle is about ± 1° as measured at the 
70% transmission points. This compares favorably with 
the calculated acceptance angle of ±1.3° at the zero 
transmission points. The two-degree acceptance angle 
gives a calculated sample volume on the order of 
1.3xlO- 2 cm 3 (13mm 3 ). 

3.2 Sample Preparation and Experimental Procedure 

Clean glassware and "particle free" water are used to 
minimize extraneous particles in the suspension. After 
each experiment, the cell is washed with detergent, 
rinsed, and then cleaned in hot nitric acid for a couple of 
hours. After rinsing the cell and cover glass 8 to 10 times 
with deionized, filtered water, the cell is dried without 
water marks by directing a particle-free flow of air over 
the external surface of the cell. Deionized water with a 



resistivity of about 18 megohm-centimeter is used for 
rinsing glassware and for sample dilution. A pleated 
membrane filter with a 0.2-^im pore size is also used on 
the water outlet to remove particulates. 

In order to minimize multiple scattering, the particle 
suspension is diluted to a concentration on the order of 
10 6 particles/cm 3 . Dilution of approximately four orders 
of magnitude is required to dilute the standard samples, 
which are 0.5% by weight. 

After positioning the optical cell, a beam stop is in- 
serted in the cell to minimize reflected light. Mea- 
surements are performed every two degrees from 20° to 
about 140°, requiring about 30 minutes per experiment. 
At each angle, two intensity measurements are made to 
check for consistency and to avoid spurious data re- 
sulting from contamination particles which occasionally 
enter the scattering volume. 

3.3 Data Analysis 

For the case of light scattering from a suspension, 
there are three unknown parameters: the number- 
average size, D n , the standard deviation of the size distri- 
bution, cr P , and the index of refraction of the poly- 
styrene sphere. We assume the size distribution to be 
normal. For the narrow size distribution of interest here, 
there is little difference between a log-normal distribu- 
tion, which others have assumed for the polystyrene 
spheres, and a normal distribution. Determinations of 
size distributions of polystyrene spheres by electron mi- 
croscopy support this assumption. In carrying out our 
data analysis, we assume the index of refraction to be 
known, and treat D a and o> as unknown parameters. In 
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subsection 3.5 we consider the sensitivity of the results 
to the assumed value of the index of refraction. 

For a distribution of particle sizes, the theoretically 
predicted scattering intensity, I(dj,D n ,n), is the single- 
particle scattering intensity, T(0i,Dj,n), averaged over 
the size distribution and the acceptance angle of the 
detector. 



I(O h D n ,cr P )= 



1 



(A -A) 2 
2o> 2 



{T(fi lt D it n) 



V(27r)<rw?. SXp 
+OJT(0 i -l,D j ,n)+O.lT(9 i +l f D s ,n)} (7) 



where m represents the number of size increments and 
the terms with 0, + 1 and 0, — 1 take into account the 
effects due to the finite acceptance angle of the detector. 
The quality of fit, Q, for the particle suspension experi- 
ment is defined by the same type of expression as eq (1), 

Q ~ | ( [E, - aI(B h D ai cr P )f, (8) 

where the index / refers to angle, and E { is the measured 
scattered intensity at 0,. The expression for a is the same 
as given by eq (2) but with 1(9,) replacing T,-. 

For the refractive index of the polystyrene spheres, 
we use the value 1.588 for bulk polystyrene at \=632.8 
nm given by Boundy and Boyer [19]. The refractive 
index of water at 23 °C was taken to be 1.3315 based on 
measurements by Tilton and Taylor [20]. The appropri- 
ate wavelength and index of refraction for a particle 
suspended in water are the ratios of the respective val- 
ues in vacuum to the refractive index of water giving 

-^=0.4752 urn 



ical experimental data and best-fit results are shown in 
figure 6. 

We also consider a second definition of quality of fit 
based on weighting the intensities by a function of the 
angle. Specifically, we consider 

/i(0,)=/(0,)sin 4 (|), (10) 

which reduces the contributions to the quality of fit 
from the large peaks at small scattering angles. As seen 
in figure 7, which is based on the same data as figure 6, 
the weighted peaks and valleys are approximately uni- 
form in amplitude throughout the angular range. The 




60° 90° 

SCATTERING ANGLE 



120° 



Figure 6-Best fit curve for the intensity of scattered light as a function 

of angle with D n = 0.893 jxm and o> =0.027 jxm. Data ( ); 

theory ( ). 



n PS L (vac) 



= 1.1925. 



The number of intervals in the size distribution, m , is 
typically taken to be 10 or 20. The size increment, A, 
from Dj to jD j+1 is chosen to be 



A = 



4o> 
m-1 



(9) 



so that the average in eq (7) covers the full size distribu- 
tion regardless of the value of o>. For a specified value 
of D n and o> , the quantity I is calculated at two degree 
intervals from 20° to 140°. The value of Q is determined 
from eqs (8) and (2). The best-fit values for D n and o> are 
determined by finding the minimum Q for a range of 
values for £> n and cr P . The range for D n is typically 
0.893-0.898 jxm in increments of 0.0005 jam, and, for a>, 
from 0.021 to 0.039 /xm in increments of 0.002 jam. Typ- 




60° 90° 

SCATTERING ANGLE 



Figure 7-Best fit curve for the intensity of scattered light weighted by 
sin 4 (0/2) as a function of angle with D„= 0.896 jxm and o>= 0.023 

jxm. The experimental data is the same as for figure 6. Data ( ); 

theory ( ). 
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values of Q as a function of D n and a> are given in table 
4 for this case (experiment A3). We consider this second 
definition of quality of fit to be superior to the first 
because all the data are significant. For our data, both 
techniques produce very nearly the same result; for ex- 
ample, D n for the first method is 0.894 compared to 0.896 
for the second method. 

3.4 Results 

One drop from one of the SRM 1690 vials, selected at 
random from 1000 vials, was used in preparing each of 
the 10 samples. Light-intensity- versus-angle data were 
taken for each sample and the best-fit parameters were 
calculated. The results of the measurements are 
presented in table 5. The average value of D n , D n , for 
experiments A2 through Jl was found to be 0.895 jam 
with a standard deviation associated with the average 
size, <x Dn , of 0.0007. The average value of the standard 
deviation of the size distribution, cr P , was found to be 
0.029 /xm. Also indicated in table 5 are four repeat mea- 



Table 4. 



The quality of fit, Q, for experiment A3 with weighting 
factor 







ar P ,ysa. 0.021 


0.023 


0.025 


0.027 


D,ju,m 










0.8945 


0.158 


0.147 


0.149 


0.165 


0.8950 


0.147 


0.140 


0.146 


0.167 


0.8955 


0.140 


0.137 


0.148 


0.172 


0.8960 


0.139 


0.140 


0.154 


0.181 


0.8965 


0.143 


0.147 


0.164 


0.192 



Table 5. Best-fit values for D n for 10 samples of SRM 1690 based on 
light scattering from a particle suspension. 











photon 




Sample" 


Aujxm 


cr 


counts/10s d 


Q 


A2 


0.8945 


0.029 


197,000 


0.276 


B 


0.8965 


0.025 


165,000 


0.287 


C 


0.8945 


0.031 


177,000 


0.321 


D 


0.8940 


0.031 


179,500 


0.272 


E 


0.8955 


0.029 


172,900 


0.288 


F 


0.8950 


0.029 


150,500 


0.224 


G 


0.8945 


0.027 


166,000 


0.246 


H 


0.8950 


0.027 


97,500 


0.166 


I 


0.8945 


0.031 


142,700 


0.231 


Jl b 


0.8955 


0.029 


137,200 


0.280 


J2 


0.8950 


0.029 


137,000 


0.221 


J3 


0.8955 


0.027 


137,500 


0.234 


J4 


0.8955 


0.029 


136,000 


0.217 


A3 C 


0.8955 


0.023 


167,500 


0.137 



"The sample concentration is approximately 7.0X 10 5 particles/cm 3 . 
b The numerals refer to repeat measurements of the same sample. 
c Sampte A3 was prepared and the measurements made six months after the 
other measurements. 
d Measured at 30 °C. 



surements for one sample with the scans performed se- 
quentially over a period of about three hours. The stan- 
dard deviation for the average size for the repeat mea- 
surements is less than 0.0003 u.m. We think that the 
larger standard deviation obtained for the 10 samples is 
a result of the sample preparation and cell positioning 
and not of differences in the vials; in any event, both 
standard deviations are small and represent a minor 
component of the overall error as discussed in sub- 
section 3.5. 

Scattering measurements were also performed with a 
He-Cd laser with X=441.6 nm, as a test for consistency. 
The appropriate values for the indices of refraction of 
water and polystyrene for this wavelength at 23 °C are 
1.3395 and 1.615, respectively. Because of the lower 
intensity stability of the He-Cd laser, the ratio of signal 
photon counts to reference counts was taken as de- 
scribed in subsection 3. 1. In this case D n was found to be 
0.896 with o>= 0.031. 

3.5 Error Analysis 

The major elements of the random error component 
of the overall uncertainty are the photon counting noise 
and the fluctuation in the number of particles in the 
scattering volume. The percentage noise arising from 
the statistical nature of photon counting goes as 1/V7V 
so that even for a low count of 10,000 the noise is only 
1%. Another 1/ViVterm is caused by the fluctuation in 
the number of particles in the scattering volume. For a 
number concentration of 7X 10 5 particles/cm 3 and a typ- 
ical scattering volume of 1.3 X 10 -2 cm 3 , the number of 
particles is 9.1 XlO 3 leading to a noise of 1.0%. As- 
suming the two fluctuations to be independent, we ob- 
tain a combined noise of 1.5%. The effect of this noise 
was simulated by the same technique as described in 
subsection 2.4 and found to change the particle size by 
0.0006 jLtm. This is comparable to the observed cr, 
0.0007, for the measurements with the 10 samples. The 
random component of the uncertainty is calculated as in 
subsection 2.4, 



R =r„_ 1 (0.025) 



Vn 



(11) 



where n refers to the number of samples. In this case the 
number of degrees of freedom is n — 1. We obtain for R 
a value of 0.0005 u-m, which is a factor of 13 less than for 
the single-particle measurements. 

The primary component of the systematic error is 
uncertainty in the refractive index of the particle with 
lesser contributions from multiple scattering, finite ac- 
ceptance angle of detector, misalignment of detector, 
and reflected light. The uncertainty in the angle mea- 
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surement and in the wavelength of light is negligible in 
comparison to these others. 

The values of the refractive index of polystyrene 
spheres obtained by various investigators were extrapo- 
lated to the He-Ne laser wavelength using the following 
dispersion equation: 



=^+T2- 



i? 

X 2 



(12) 



The constants A and B are obtained from the measured 
value of n at one wavelength together with Starkie's [21] 
value of 31.0 for the dispersion ratio defined as 



^5893 — 1 
n4861 — ^6563 



■=31. 



(13) 



The resulting values given in table 6 are close to the bulk 
value of 1.588 but, based on these data, we cannot rule 
out the possibility of the refractive index of the particle 
being slightly different from that of the bulk. The range 
of values obtained for liquid suspensions of the particles 

Table 6. Refractive index of polystyrene. 





Investigator 


n 


\,jxm 


n(0.6328) a 


Boundy and Boyer 








bulk value 






1.588 


Heller and Pugh 








0.132-0.295 ju,m 


1.602±0.008 


0.5461 


1.594 


Smart and Willis 








0.132-0.295 jam 


1.581 


0.589 


1.577 


Marx and Mulholland 








D = 0.460 jam b 


1.603 ±0.029 


0.4416 


1.579 


D= 0.600 jxm 


1.599±0.013 


0.4416 


1.575 


D=0.914 jum 


1.612±0.004 


0.4416 


1.588 


D = 1.101 u.m 


1.619±0.002 


0.4416 


1.595 


McRae 








D =0.600 jum b 


1.63±0.019 


0.5145 


1.618 


D= 0.794 fun 


1.58±0.008 


0.5145 


1.569 


D= 1.011 jum 


1.57±0.012 


0.5145 


1.559 


Davis and Ravindran 








D=11.9 and 5.7 fim 






1.575 


Bottiger, Voss, 








and Fry 








D=1.112 /xm 


1.615 


0.4416 


1.590 


Phillips, Wyatt, 








and Berkman 








D= 1.099 ftm b 


1.59±0.001 


0.5145 


1.578 



The refractive index has been extrapolated to X =0.6328 using the equation 
n=A+B/\ 2 . 

"The values of D are the values obtained by Dow Chemical Co. using trans- 
mission electron microscopy. 



extends from 1.577 to 1.594 and is based, respectively, 
on transmitted light intensity measurements by Smart 
and Willis [22] and interferometric measurements by 
Heller and Pugh [23]. In both cases the data analysis is 
the most straightforward for small particle sizes, and 
only the small particle size results are given in table 6. 

All the other values in the table correspond to light- 
scattering measurements of single particles levitated in 
air. We consider Marx and Mulholland's value [14] of 
1.588 for the nominal 0.914 u.m particles to be the best of 
the single-particle measurements in regard to SRM 
1690, since the measurements were made on the SRM 
particles and great care was taken with regard to angle 
calibration and data analysis. The low value obtained by 
McRae [11] might have resulted from an error in the 
calibration of the angle encoder. We observed an appar- 
ent size-dependence of the refractive index similar to 
that reported by McRae when our angle calibration was 
off by about 3°. The measurements of Davis and Rav- 
indran [12] were done on much larger polystyrene 
spheres (5.7 and 11.9 jim) and the quality of the data 
appears to be lower than that of the other studies. The 
poorer data may have resulted from a slight asphericity 
or inhomogeneity in the particles. We take as our range 
of uncertainty in refractive index 1.577 to 1.595, which 
encompasses all the values in table 6 except the values of 
Davis and Ravindran, McRae, and Marx and Mul- 
holland's value for the 0.600- urn particle size which had 
a large uncertainty. 

To estimate the effect of this uncertainty on D n , we 
repeated the data analysis for experiment A3 using the 
two extreme values of n. The value 1.577 leads to a 
slightly increased value of D n , 0.899, while the larger 
value of 1.595 leads to a slightly reduced mean diameter 
of 0.893 jam. So the contribution to the uncertainty from 
this effect, S„, equals ±0.003. 

The effect of multiple scattering on the results was 
estimated by performing measurements at two concen- 
trations. According to the criterion established by 
Napper and Ottewill [24], multiple-scattering effects be- 
come insignificant for the ratios of the average distance 
between particles, /, and the particle diameter, D, 
greater than 100. For concentrations of 7.0 XlO 5 and 
3.5 X 10 6 particles/cm 3 , we find that l/D is 125 and 70 
respectively. From analyzing the light-scattering data 
for these two cases, we obtain 0.894 \im at the lower 
concentration and 0.890 jam at the higher concentration. 
We estimate 8 m as +0.001 fxm allowing for a slight 
residual effect even though the criterion given above is 
satisfied for the 10 samples analyzed. 

The effect of the finite acceptance angle of the pho- 
tometer was determined by computer simulation in 
which a slit width integration was performed over the 
theoretically-predicted intensity pattern for D n = 0.900 
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and o> =0.010 jam. The "data" were then analyzed by 
the same procedure with the intensity weighted by 
sin 4 (0/2), and the resulting values of D n and cr P were 
0.900 and 0.011 jam respectively. So we find 8 e to be at 
most 0.0005 jam. If the intensities were analyzed directly 
without the angular weighting function, the values of D n 
and cr P , 0.899 and 0.015, changed more as a result of the 
finite acceptance angle than the values obtained using 
the intensity weighted by the angular function. We also 
used computer simulation to estimate the error associ- 
ated with a slight misalignment at zero angle, about 
0.07°, and found 8 D to be 0.0004 jam. 

Without a beam stop, on the order of 4% of the laser 
beam is reflected from the glass-air interface back into 
the cell. This would lead to some fraction of the light 
reaching the detector having been scattered through the 
supplement of the scattering angle. We obtain an in- 
crease of 0.002 jam size without the beam stop and esti- 
mate an error, 8 R , of at most +0.001 jam due to residual 
reflection. 

As discussed in section 5, about 1.5% of the particles 
is in the form of agglomerated doublets. Bottiger et al. 
[13], having measured all 16 components of the Mueller 
matrix for a doublet of 1.112 jam polystyrene spheres, 
find that the peaks and valleys in the scattering pattern 
appear at the same angles as for the single sphere but 
that the amplitudes are smaller. We have no quantitative 
estimate of this agglomeration error, 8 A , but we intu- 
itively expect it to be less than 0.001 jxm. 

Estimating the total uncertainty in the same manner as 
in subsection 2.4, we obtain an expression for U T shown 
below together with estimates of the different terms, 

C/ T =i? + |8 n | + |S m | + |8 fl | + |8 D | + |0 R | + |8 A |. 
0.0005 0.003 0.001 0.0005 0.0004 0.001 0.001 



Using these estimates, we obtain C/ T = 0.007 jam with a 
number-average size, D n , of 0.895 jam. The uncertainty 
is lower than for the single-particle measurement be- 
cause of the much smaller random component of error, 
0.0005 jam compared to 0.007 jam. 

While the average size can be determined quite accu- 
rately by scattering from a suspension, this is not the 
case for the standard deviation of the size distribution, 
<r P . Our best estimate of a> based on light scattering is 
0.029 jam while, as discussed in section 5, the value 
obtained by electron microscopy is about 0.0095 jam. 
We suspect that the major cause of the discrepancy is 
the contribution of scattered light from off-size particles 
and agglomerated doublets. The off-size particles are 
excluded from the determination of cr P by electron mi- 
croscopy and the doublets are treated as pairs of individ- 
ual spheres. 



4. Optical Array Sizing 

Array sizing by optical microscopy provides an inde- 
pendent measure of D n since the technique is based on 
geometrical optics rather than diffraction theory. This is 
one of the oldest techniques for the measurement of 
average size of micrometer-sized particles. Perrin [25] 
used this technique in his study of Brownian motion, 
from which Avogadro's number was first determined. 
Its ability to produce accurate values for the average 
diameter results from a combination of a highly accurate 
length standard for optical microscopy and the fact that 
monosize polystyrene spheres can form close-packed, 
two-dimensional hexagonal arrays. The microscopic im- 
age of a row of micrometer-size spheres consists of a 
string of circular images blurred by diffraction as indi- 
cated in figure 8; however, the center of the sphere can 
be pinpointed with considerable accuracy. In fact, when 
rows of, say, 10 spheres are used, the potential exists for 
finding the average size of 1 jam spheres to a resolution 
of about 1% of the sphere diameter. Our methods for 
preparing the arrays and performing the measurements 
are similar to those described by Kubitschek [26]. 

4.1 Formation of Two-Dimensional Arrays 

To minimize non-volatile residue, the particles were 
"washed" with 18-megohm-centimeter deionized water 
filtered with a 0.2 jam pore size pleated membrane filter. 
The washing consisted of pouring off the supernatant 
after allowing the particles to settle and then adding 
filtered, deionized water. The particle concentration 
was about 0.5% by mass. A drop of the particle sus- 
pension was placed at one end of a microscope slide and 
smeared out by means of a second slide held such that 
one of its long edges would wipe slowly over the first 
slide. The evaporating water sweeps the particles in- 
ward though it appears that surface tension is re- 
sponsible for the small scale motion leading to the or- 
dered array. In the early phase of evaporation, the 
particles tend to congregate in assemblies similar in 
shape to ice floes in a freezing ice field. As the evapo- 
ration continues, the spheres arrange themselves into 
well-ordered hexagonal arrays, with rows 10-50 spheres 
long. This ordering occurrs too rapidly to be observed 
by eye. 

The ordered arrays extending perhaps 7 spheres in 
one direction and 20 in another are separated by micro- 
cracks as seen in figure 9. The cracks occur during the 
drying when an array field is pulled apart perhaps to 
relieve a stress on the array. Another type of crack, also 
illustrated in figure 9, is much shorter and is V-shaped. 
It seems to indicate intensive drying within an array 
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Figure 8-Hexagonal array of nomi- 
nal one-jnm polystyrene spheres 
(Standard Reference Material 
1690). 



field. Close inspection shows that such cracks can be 
detected when they are about 0.05 u.m wide, and they 
are easily visible when the width exceeds 0. 1 jam. 

4.2 Measurement of Arrays 

The arrays were measured by scaled micro- 
photography rather than by microscopy with a gradu- 
ated eyepiece. Microphotography has the advantage of 
providing a permanent record from which array defects 
such as microcracks can be accounted for. The image 
distortion in the optical train was measured to be less 
than 0.2%, and scale changes from repeated focusing 
were 0.1-0.2%. The microscope was a Zeiss Universal 1 
with Polaroid attachment and with a 63 X, 0.90 N.A. 
objective. 



The micrographs, with a typical scale of 1.6 mm per 
u.m, were measured to 0.05-0. 1 mm. The magnification 
calibration and the image distortion were measured us- 
ing a 19.990 ±0.005 ju-m linewidth spacing on SRM 474 
Line Width Standard from NBS. 

The arrays were prepared on six microscope slides. 
One to three arrays were measured on each slide for a 
total of 12 arrays. For each array, 15 row lengths were 
measured. The number-average diameter for a single 
row is determined by dividing the length of the row, L , 
by the number of spheres in the row minus one, N—\. 



A= 



N-\ 



(14) 



The —1 results from using the center-to-center row 
length rather than the edge-to-edge row length. 
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Figure 9-The same field of view is slightly defocussed for the photograph on the right to allow better visualization of the microcracks. 



In our measurements the row length varied from 10 to 
30 spheres. 

4.3 Results 

The results for two arrays are shown in table 7. Mea- 
surements along each of the three lattice directions are 
included. In some cases the row lengths in different 

Table 7. Measured row lengths for two arrays. 

Orientation 
Row Length Readings Row Average Diameter Average 
(mm) (jLim) (jam) 



Slide 








No. 


0° 


120° 


240° 


3 


28.6 


29.0 


29.15 




28.8 


29.0 


29.15 




28.9 


28.0 


29.2 




28.85 


28.95 


29.1 




28.9 


28.9 


29.1 


6 


14.55 


14.55 


14.7 




14.5 


14.55 


14.7 




14.5 


14.5 


14.7 




14.5 


14.4 


14.7 




14.5 


14.55 


14.6 



0° 



120° 240° 



0.891 0.901 0.896 



0.896 



0.898 0.898 0.908 



0.901 



directions differ by as much as 0.2 mm for an 11-sphere 
array leading to differences in the average diameter of 
0.01 jam. The reason for this anisotropy is not clear. 
Kubitschek [26] attributes this effect to undetected mi- 
crocracks. If the cause were a single microcrack, then it 
would be on the order of 0. 1 to 0.2 mm and would thus 
have been detected. For those arrays located at the pe- 
riphery of the ordered region, we found that in general 
the average diameters calculated from rows parallel to 
the edge were somewhat larger than values calculated 
from rows in the other two directions, suggesting loose- 
ness in that direction. The drying speed of the array 
also seemed to have some influence. 

The row-average diameter, D„ is obtained as a simple 
average of the values D T (J) for the R rows. 



2A(0 
Ux R 



(15) 



The standard deviation in D t ,cr Dt is given by 



2(A(0-A) 2 



<r D =^ 



R-l 



(16) 



A summary of the values of D t and cr Dx for each of 12 
arrays, which includes 240 rows, is given in table 8. The 
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Table 8. Array sizing results for SRM 1690. 



Slide 


Array 


Row 


No. of 


A 


°"D r 


Anisotropy 


No. 


No. 


Length 


Rows 


[im 


firn 


fxm 


1 


1 


20 


15 


0.905 


0.0027 


0.001 


1 


1 


20 


15 


0.909 


0.0027 


0.006 


1 


1 


20 


15 


0.906 


0.0037 


0.003 


1 


1 


20 


15 


0.905 


0.0018 


0.002 


1 


1 


20 


15 


0.904 


0.0028 


0.002 


1 


2 


15 


15 


0.895 


0.0027 


0.005 


1 


3 


15,30 


15 


0.896 


0.0027 


0.005 


2 


4 


20 


15 


0.896 


0.0022 


0.002 


3 


5 


20 


15 


0.898 


0.0049 


0.011 


4 


6 


13 


15 


0.886 


0.0057 


0.012 


4 


7 


13 


15 


0.888 


0.0038 


0.009 


4 


8 


13 


15 


0.889 


0.0039 


0.007 


5 


9 


20 


15 


0.882 


0.0029 


0.004 


6 


10 


20 


15 


0.898 


0.0023 


0.003 


6 


11 


10 


15 


0.903 


0.0044 


0.010 


6 


12 


10 


15 


0.902 


0.0039 


0.010 






Average: £>„= 


=0.898 


0.0026 


0.006 








er Dn : 


=0.008 







average of all the values of D T in table 8, defined as D n , 
is 0.898 with a standard deviation among the row aver- 
ages, <r Dn , of 0.008 jam. The anisotropy, defined as the 
range in the row-averaged diameters in the three lattice 
directions of an array, is found to have an average value 
of 0.006 jam. This is comparable in magnitude to cr Dn . 
It is also possible to determine the standard deviation 
of the size distribution, cr P , from cr Dn . Assuming the 
spheres to be randomly distributed along a row in terms 
of particle size, it can be shown that the cr P is related to 
cr Dr by the expression, 



o>=ViV— 1 cr Dr , 



(17) 



where N is the number of spheres in a row. 

In arriving at this result, one makes use of the fact that 
the variance of the sum of random variables (the row 
length) is the sum of the variances of the individual 
random variables (sphere diameter). The resulting 
value, 0.013 jam, has a statistical uncertainty of about 
±0.003 jam and is somewhat larger than the value ob- 
tained by electron microscopy, 0.0095 firn. 

4.4 Error Analysis 

While the experimental technique of optical array siz- 
ing is relatively simple compared to light-scattering 
measurements, there are more potential sources of un- 
certainty, both random and systematic. The random 
component of the total uncertainty is related to the stan- 
dard deviation of the average diameters for the arrays 



measured, <r Dn . We estimate the magnitude of the vari- 
ous components of the random error including image 
distortion, readout of microphotographs, local motion 
in the photographic material, and anisotropy. These are 
combined and compared with cr Dn . There are also a 
number of systematic errors arising from the finite width 
of the size distribution, microcracks, solute impurities, 
and shape distortion. 

Image distortion results in the length scale changing 
over the photographic field of view. Because of axial 
symmetry of the optical train, it needs to be measured in 
relation to one plane containing the optical axis. The 
Standard Reference Material 474 Line Width Standard 
was used as a two-line resolution target. 

A series of 150 measurements was made, covering 
five photographs, each with 10 areas and each area be- 
ing measured three times. The result was a scale factor 
of 0.6 19 ±0.001 jam/mm (— at center, + at edge). Using 
one scale value for the entire field of view introduces an 
image distortion error, cr u no greater than 0.002 jam. 

Row lengths were measured to 0.05-0.1 mm 
(0.03-0.06 urn) using a film scale. The rows contained 
10-20 spheres; thus the error contribution to the average 
diameter value (1.0 firn nominal) was 0.003 to 0.006 jam. 
In practice the film readout error, o- 2 , was held below 
0.005 jam. 

It is known that slight motion can occur in photo- 
graphic materials primarily as a result of uneven drying. 
In the Polaroid material used here, there was no de- 
tectable motion at the 0.05 mm level. The error con- 
tribution, o- 3 , is thus less than 0.003 jam for 20-sphere 
rows. 

The anisotropy in the array will also manifest itself as 
a component of the random error, o- 4 . We estimate this 
effect as half the range in the average anisotropy, 
0.003 jam. 

We obtain a combined estimate of all the sigmas by 
using the formula 



cr t 2 = 2 oV 



(18) 



In using this standard formula, we are assuming the 
errors to be independent. The resulting value of cr t , 
0.007 jam, is comparable to or Da so our error budget 
accounts for the observed value of the random error. 
The random component of the overall uncertainty based 
on 16 sets of measurements is found to be 0.0044 jam 
based on the formula 



R =f 15 (0.025) 



<r Da 
Vl6 



(19) 



It is thought that the mean distance between particle 
centers in monodisperse arrays, such as those shown in 
figure 8, is greater than the mean particle diameter be- 
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cause many smaller particles fail to be in contact with all 
of their neighbors. Thus, gaps in the array add to the 
length of a line of particles, causing an overestimate of 
their true mean diameter. For hexagonal arrays formed 
from machined washers of known size, Kubitschek [27] 
found that the overestimate AD for particles with stan- 
dard deviation o> was given by 



AZ>=0.46o- f 



(20) 



For our case, cr P obtained from electron microscopy is 
0.0095 jam, so the resulting systematic error, 8 W , is 
—0.0044 jam. We have chosen to treat this effect as a 
systematic error rather than correcting for it, because 
there has been no direct verification of eq (20) for 
micrometer-size spheres. 

A second source of systematic error arises from mi- 
crocracks. When a row is measured across one or two 
undetected cracks, one can expect an error contribution 
to the average diameter of about 0.01 jxm for rows 10 
spheres long, and less for longer rows. We take as our 
error estimate, 8 C , a value of —0.005 jam. 

The presence of impurities in the polystyrene sus- 
pensions is another source of error. Assuming that the 
emulsifier and inorganic impurities form a uniform coat- 
ing around the spheres when the water evaporates, one 
finds the following expression for the thickness, t, of the 
coating 



t = 



C\D 
6C 



(21) 



where C is the volume concentration of the polystyrene 
and Q the concentration of the non-volatile impurities. 
For our case the non-volatile impurity concentration 
was about two orders of magnitude lower than the vol- 
ume concentration of the polystyrene spheres, and the 
maximum resulting error is about 0.003 jam. 

An attempt was made to detect this effect by using 
one sample of the stock solution and a different one in 
which the original solute impurity content had been 
reduced by an order of magnitude by settling, decanting, 
and addition of filtered, deionized water. The measured 
diameter for the second sample was slightly less, by 
0.004±0.003 jam. Therefore, all reported measurements 
were performed on samples with reduced impurity con- 
tent. We estimate the error resulting from residual non- 
volatile impurity of washed particles, 8 b as 0.001 jam. 

Two other potential systematic errors include the ef- 
fect of unequal-size spheres on a flat substrate and the 
effect of curved rows. If we have a row of spheres of 
unequal diameters, the centers of the spheres will not be 
on a straight line and the row length will be less than the 



sum of the diameters. If the size distribution is not too 
broad (<r P /D <0.03), the error in the calculated average 
diameter, cr P 2 /2D 2 , is small (less than 0.001 jam) com- 
pared to other errors. 

Occasionally arrays will exhibit curved or wavy 
rows. For typical curvatures, replacing the arc by the 
chord introduces only a small error (on the order of 
0.001 jam) in the measured diameter. However, we 
found that curved rows are sometimes as much as 1% 
longer than the corresponding straight rows found in 
the same array. This suggests the presence of slack in the 
curved rows, and therefore such rows were omitted 
from our measurements. 

It is possible that spheres in contact flatten slightly. 
We made an attempt to measure the average size of 
isolated spheres and arrays of spheres with a trans- 
mission electron microscope, but we were not able to 
prepare arrays on TEM grids. There were instances of 
particles in contact as in figure 10, and from measuring 
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Figure 10-Electron micrograph of nominal one-jxm polystyrene 
spheres (SRM 1690) at a magnification of 10,000 X. The spacing for 
the diffraction grating lines is 0.46 jum. 
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such particles we found that the flattening effect is less 
than 0.005 jam. More quantitative measurements are 
necessary to refine 8 F . 

As in the error analyses in the previous sections, we 
obtain the combined uncertainty from the sum of the 
random error, R , and the absolute values of the system- 
atic errors, 

C/ T = J R + |8w| + |8 c | + |5 I | + |S F | 

We obtain a value of 0.020 for U T with the overall 
average value for D n of 0.898 juim. 



5. Measurements of Size 
Distribution/Concentration 

We have used transmission electron microscopy, 
flow-through electrical-sensing zone counter mea- 
surements, and optical microscopy to obtain more de- 
tailed information about the size distribution. Using 
these techniques, we measure a>, the standard deviation 
of the size distribution, the fraction of off-size particles, 
and the fraction of agglomerated doublets. Our only 
direct measure of number concentration was obtained 
with electrical sensing zone counter measurements. 

5.1 Transmission Electron Microscopy 

The sample was prepared for electron microscopy by 
evaporating a small drop of diluted suspension of the 
spheres on a diffraction grating replica mounted on a 
TEM grid. Three grids were prepared from one SRM 
1690 vial. The grids were coated with approximately 20 
nm of carbon and then examined with a JEOL 200CX 
electron microscope 1 at an accelerating potential of 100 
KV and at a nominal 10,000 X magnification. The par- 
ticle size was measured directly from the negative using 
a 7X magnifier with an accurate millimeter scale reti- 
cule. A series of concentric circles on the reticule en- 
abled quick location of the sphere diameter. A typical 
micrograph is shown in figure 10. 

As mentioned in the introduction, we had hoped to 
make an accurate determination of the average particle 
size from the TEM measurements. However, we found 
significant variability in the commercially available gra- 
ting replicas in regard to the line spacing even among 
those from the same supplier. We also found variations 
as much as 3% in the line spacing in different regions of 
the same grating. So from our assessment, transmission 
electron microscopy is not competitive with the light 
scattering and array sizing measurements for accurate 
particle size measurement without a suitable calibration 
artifact. 



However, the high resolution of the TEM makes it 
suitable for the measurement of the standard deviation 
of the narrow size distribution. The size distribution 
typically consists of a very sharp central peak plus 
broadly distributed off-size particles. The inclusion of 
the off-size particles in the data analysis will greatly 
broaden the apparent value of o>. In one case cr P 
changed from 0.031 to 0.011 as a result of removing 
three off-size particles from a sample of 102 particles. It 
is, therefore, important to have a systematic procedure 
for rejecting off-size particles. We used a discordancy 
test [28] based on the sample kurtosis as the test statistic, 



sample kurtosis =- 



*■ ' i'= 1 

CJCA-A) 2 ) 2 ' 



(22) 



If the sample kurtosis exceeds a value of about 3.77 for 
a 100-sphere sample, then one or more spheres are off- 
size at the 5% level of discordancy. The spheres with 
diameters farthest from the average size are eliminated 
consecutively until the sample kurtosis reaches the ap- 
propriate value. The particles at the edge of the field of 
view were eliminated because of evidence of mag- 
nification distortion in this region. 

Using a common magnification for all the micro- 
graphs resulted in a smaller value of o> than the value 
obtained using a magnification for each photograph de- 
rived from the line spacing on the diffraction grating 
replica. The smaller value of cr P is the better value be- 
cause experimental errors only tend to increase cr P . The 
constancy of the magnification was to be expected, since 
all the micrographs were taken at one magnification 
setting over a period of a few hours. We obtained values 
of 0.0083 jam and 0.0106 /xm for o> based on samples of 
59 and 99 particles. We take the mean value of 0.0095 as 
our estimate of o>. 

The major source of uncertainty in o> is statistical in 
origin. The sample variance has a x 2 distribution, and for 
a sample size of 100 there is an uncertainty of about 
±13% at the 95% confidence level. The size distribu- 
tion is represented as a histogram in figure 11. 

Other information obtained from TEM includes the 
asphericity of the particles and the fraction of "new" 
particles. The asphericity, A„ is defined as 



"( a 



-Dn 



D 



X 100%. 



(23) 



The diameter measurements were made every 45° for 
particles located near the center of the field of view. The 
asphericity is small, 0.6% ±0.3%, with large uncertainty 
because of the resolution limit of the TEM. 

We found a total of eight particles with diameters less 
than 0.5 am out of some 1350 particles (0.6%). As ex- 
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Figure 11-Size distribution of nominal one-urn polystyrene spheres 
(with off-size particles excluded) obtained by transmission electron 
microscopy. The dashed curve is the normal distribution with D n 
(0.932 fim) and o> (0.0 106 fim) obtained from the experimental data. 

plained in the Appendix, it is possible to nucleate a new 
population of particles in the process of using seed par- 
ticles to produce larger particles. We believe that this 
nucleation process is responsible for the small off-size 
spheres. 

5.2 Electrical Sensing Zone-Counter Measurements 

We used the Coulter ZB Counter with Coulter Chan- 
nelyzer Pulse Sizer 1 to measure the fraction of doublets, 
the fraction of large off-size particles, and the number 
concentration, and to estimate the variability of D n in the 
SRM 1690 samples. To apply this technique, the par- 
ticles must be prepared in an electrolyte solution at a 
concentration of about 1% by mass. The technique con- 
sists of measuring the electrical current due to the elec- 
trolyte flowing through an orifice; in our case, 30- and 
19-jam diameter orifices were used. As a particle flows 
through the orifice, the current is reduced by an amount 
approximately proportional to the volume of the par- 
ticle [29]. 

For our application, we used nominal 0.9- and 2-jam 
size polystyrene spheres for a particle size calibration of 
the instrument, and then used the Coulter Counter to 
measure the concentration and size distribution. The 
Coulter Counter is not an absolute sizing instrument. 

It was found that if the polystyrene spheres were 
directly diluted with electrolyte, then about 10% ag- 
glomerated doublets resulted from the combined effects 
of high particle concentration and the reduction in elec- 
trostatic repulsion between the particles. The agglomer- 



ated doublets were minimized by first diluting with fil- 
tered, deionized water (dilution factor 2.5 XlO 3 ) and 
then diluting by another factor of 10 with the electrolyte 
solution. The percentage of doublets was found to be 
1.5 ±0.4% of the total population for fresh samples, 
compared to about 1 % obtained by visual microscopy of 
the suspended particles diluted without electrolyte. 
Based on the limited sample size, there does not appear 
to be any significant difference in the fraction of dou- 
blets between the sample diluted without electrolyte 
and the fresh samples diluted with electrolyte at the last 
step. 

By changing the electronic settings of the Coulter 
Counter, it was possible to estimate the percentage of 
large single particles with diameters in the range of 2— 6 
fxm. The observed count was on the order of 200 out of 
a total population of about 235,000 particles, which cor- 
responds to slightly less than 0.1%. 

Based on the measured number-concentration of two 
diluted samples for particles in the 0.6 to 1.4 fxm size 
range, we obtain number concentrations of the un- 
diluted samples of 1.18 and 1.45 XlO 10 particles/cm 3 . 
These values agree well with a calculated value of 
1.25 XlO 10 particles/cm 3 based on a concentration of 
0.5% by mass of the polystyrene spheres, a particle den- 
sity of 1.052 g/cm 3 , and a particle size of 0.9 jam. 

Diluted samples were prepared from each of 10 SRM 
vials and two portions of each sample were sized with 
the Coulter Counter using the 30-jLtm-diameter orifice. 
A typical size distribution based on sizing about 85,000 
particles is given in figure 12. A complete listing of the 
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Figure 12-Size distribution of SRM 1690 obtained by a Coulter Coun- 
ter for aperture current, amplification, and lower threshold settings 
of 1/3, 1/2, and 100 respectively. The instrument output is linear 
with respect to particle volume and consequently appears stretched 
when plotted versus particle diameter. 
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average particle size, D n , and the standard deviation of 
the central peak are given for each sample in table 9. The 
standard deviation associated with 10 independent de- 
terminations of D n has a value of 0.0027 jam. This value 
is approximately three times greater than the value ob- 
tained by light scattering and is due to instrumental 
noise rather than sample-to-sample variability. Even 
when measurements are made on the same sample such 
as Fl to F4, variations in D n as large as 0.005 jam are 
observed. 

Table 9. D„ and cr P for Coulter Counter measurements of SRM 1690. 





Portion 1 


[ 


] 


Portion 2 


Sample 


D„, fun 


crp, /xm 


A, /xm 


crp, jum 


A 


0.9226 


0.0370 


0.9239 


0.0336 


B 


0.9225 


0.0380 


0.9240 


0.0344 


C 


0.9266 


0.0351 


0.9286 


0.0349 


D 


0.9226 


0.0336 


0.9240 


0.0339 


E 


0.9264 


0.0351 


0.9251 


0.0357 


Fl a 


0.9308 


0.0392 


0.9302 


0.0332 


F2 






0.9251 


0.0329 


F3 






0.9246 


0.0326 


F4 






0.9245 


0.0325 


G 


0.9245 


0.0328 


0.9247 


0.0356 


H 


0.9231 


0.0333 


0.9222 


0.0338 


I 


0.9223 


0.0338 


0.9236 


0.0340 


J 


0.9226 


0.0323 


0.9228 


0.0311 



The numerals refer to repeat measurements of the same sample. 

The average value of the standard deviation of the 
size distribution, 0.033 /Am, for the central peak based on 
Coulter Counter measurements is more than a factor of 
three greater than the value obtained by electron mi- 
croscopy. The primary cause of this broadening is the 
variation in the electrical pulse resulting from different 
particle trajectories through the orifice. The amplitude- 
time history of the electrical pulse for a particle moving 
near the aperture wall is different from that of a particle 
moving near the centerline. Thorn [30] demonstrated 
this effect on a model system with plastic washers and 
Spielman and Goren [31] demonstrated that the resolu- 
tion of the Coulter Counter could be improved by hy- 
drodynamically focusing the particles through the cen- 
ter of the orifice. The Coulter Counter we used did not 
have hydrodynamic focusing, but it did have a pulse 
analyzer for rejecting counts in part due to particles 
with atypical trajectories. 

6. Summary 

The values of D n and the associated uncertainty for 
three techniques are: 0.895 ±0.007 /xm for light scatter- 
ing from particles suspended in water, 0.900+0.012 jam 
for single-particle light-scattering measurements, and 
0.898±0.020 jam for optical array sizing. For each tech- 



nique there is a consistency test in addition to a quan- 
titative error analysis. For the single-particle light- 
scattering measurements, the particle diameter obtained 
for the incident light polarized in the vertical direction 
(0.903 /xm) agreed well with the diameter obtained for 
light polarized in the horizontal direction (0.898 jam). 
For light scattering from a particle suspension, the par- 
ticle diameter obtained with two different laser sources, 
He-Ne (0.633 jam) and He-Cd (0.442 jam) differed by 
only 0.001 jam. In the case of array sizing, six indepen- 
dent samples were prepared and the variation in the 
average particle diameter, 0.882 to 0.903 jam, was less 
than the overall uncertainty. 

We use 0.895 jam as the certified value because of the 
lower uncertainty associated with light scattering from 
a suspension. Another advantage of this value is that the 
measurements characterize the particles in the same 
form as provided in the SRM 1690, this is, in suspension 
rather than as a dried sample where there is the possi- 
bility of a coating of non-volatile impurities on the sur- 
face. It is reassuring that the variability among the three 
values of D n is less than the uncertainty associated with 
any one technique. We expect that our estimates of un- 
certainty are, if anything, conservative. The value ob- 
tained by Dow Chemical Co. for D n by TEM, 0.914 jam, 
is about 2% greater than our certified value. 

In table 10 we summarize all the information regard- 
ing the size distribution obtained by our measurements. 
Only for D a and the index of refraction have detailed 
error analyses been made. The different features of the 
global size distribution extending from 0. 1 to 6 jam are 
indicated qualitatively in figure 13. 

A wide range of values is found for the standard devi- 
ation of the size distribution, cr P , extending from 

Table 10. Summary of results for SRM 1690. 

A,u.m 0.895 ±0.007 a light scat, from suspension (LSS) 
0.900±0.012 single part, scattering (SPS) 
0.898±0.020 opt. array sizing (OAS) 

crr.nm 0.0095±0.002 TEM, 0.010 SPS, 0.013 OAS, 0.029 LSS, 
0.033 Coulter Counter 



% doublets 


1.5±0.4 


Coulter Counter 


% large off-size 


0.1 


Coulter Counter 


particles 






% "new" particles 


0.6±0.2 


TEM 


% asphericity 


0.6±0.3 


TEM 


number concentration 


1.3+0.2 X10 10 part/cm 3 


Coulter 


index of refraction 


1.612±0.009 


single particle 
scattering 


sample variability 


0.0005 /xm 


light scat, from 
suspension 




0.0009 /xm 


Coulter Counter 



a The first value in each case is the best estimate. 



22 



Figure 13-Global size distribution. 
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0.0095 )Ltm to 0.033 jLtm. This is because each instrument 
has a sigma associated with repeat measurements of the 
same objects, cr r , and measured values of cr P given in 
table 10 represent a combined effect of the cr of the true 
size distribution and cr t . The single-particle scattering 
measurement has the smallest cr r , 0.002 jam, but the sta- 
tistical uncertainty due to the small number of particles 
measured (eight) limits the accuracy of ar P . The best 
estimate of cr comes from TEM with a slightly larger 
cr„ about 0.003-0.004 /xm, compared to the single par- 
ticle measurements, but this is offset by the improved 
statistics from a much larger sample size (100). The large 
values of cr T for the other instruments result from the 
dependence of the electrical pulse on the particle tra- 
jectory through the orifice for the Coulter Counter, the 
presence of doublets and off-size particles for light- 
scattering measurements of a particle suspension, and 
array nonuniformities for the case of optical array siz- 
ing. 

Each measurement method requires a slightly differ- 
ent sample preparation. A necessity in almost all cases is 
particle free water for dilution. We found 18 megohm- 
cm deionized water passed through a 0.2 /xm pore size 
filter to be adequate for our measurements. Clean glass- 
ware is a necessity for diluting the samples and a clean 
optical cell free of water marks is required to obtain 
good light scattering data. As described in section 3, we 
used both detergent and hot, concentrated nitric acid for 
cleaning our glassware. In the case of diluting a sample, 
especially if an electrolyte is used, there is always the 
possibility of promoting agglomeration through the re- 
duction in the electrostatic repulsion between the par- 
ticles. We found that even for a 1% by weight solution 
of sodium chloride little agglomeration occurred in an 
hour for a particle concentration of about 10 6 part/cm 3 
if the electrolyte were introduced only at the final dilu- 
tion stage as described in subsection 5.2. 



7. Concluding Remarks 

Our measurements have resulted in a particle size 
standard with an uncertainty in particle diameter of 
about 0.8%. Modest improvements in sizing accuracy, 
say a twofold reduction in uncertainty, may be possible 
with the existing techniques. About half the total uncer- 
tainty in the scattering measurements for suspensions is 
associated with the uncertainty in the refractive index. 
The refractive index of liquids is known io an uncer- 
tainty of a few parts in 10 5 , while the uncertainty in the 
refractive index of the polystyrene spheres is a few parts 
in 10 3 . It is possible that at least a tenfold reduction in the 
uncertainty in the refractive index of the spheres could 
be achieved by index matching the spheres with a liquid. 
A second improvement would be to remove the ag- 
glomerated doublets by sedimentatioi before per- 
forming the scattering measurements. 

An accurate value of the refractive index of the poly- 
styrene spheres would also improve the sizing capabili- 
ties of single particle scattering messurements since 
there would be only one unknown parameter, the par- 
ticle diameter, to be determined from the scattering 
measurements. We have found that tie sizing resolution 
improves by a factor of two to thrse if the refractive 
index is fixed when performing the data analysis. The 
single greatest source of uncertaiity arises from the 
small number of particles sized (ei,'ht); this component 
of the error could be cut in half if 4* particles were sized. 

Half of the uncertainty for optical array sizing results 
from smaller particles not toucling neighboring par- 
ticles and from slight flattening of the spheres on con- 
tact. Hartman [32] has demonstflted that it is possible to 
make square arrays of particles, where, unlike in- 
hexagonal arrays, every partite is in contact with its 
neighbors. A systematic analysis of such arrays is possi- 
ble but would be time-consuming because of the rela- 
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tively rare occurrence of square arrays. A quantitative 
analysis of flattening could be made by transmission 
electron microscopy if arrays could be prepared on 
TEM grids. 

Major improvements in sizing accuracy, say a tenfold 
reduction in uncertainty, will require a more detailed 
characterization of the particle structure. Specific struc- 
tural information of interest includes particle shape, es- 
pecially in regard to deviation from sphericity, particle 
inhomogeneity in regard to a variable index of refrac- 
tion within the particle, and the surface structure both 
for particles in suspension where an electrical double 
layer exists and in aerosol form where nonvolatile im- 
purities coat the surface. The pertinent structural infor- 
mation must also be included in the light scattering the- 
ory for determining particle size from scattering 
intensity as a function of angle. Of course, if particles 
could be made with a high degree of sphericity and low 
inhomogeneity within the particle, the analysis would 
be simplified. 

Ultimately it is planned to develop particle size stan- 
dards over the diameter range of 0.1 to 100 jam. For 
such a wide range, a variety of techniques will be used. 
One promising instrument for larger particle sizes in- 
volves optical interferometry coupled to an electron 
microscope. One of us (G.G. Hembree) has developed 
such an instrument with a piezoelectric stage, the dis- 
placement of which is monitored interferometrically 
with an uncertainty on the order of 0.02 jam. Another 
technique for larger particles consists of the mea- 
surement of resonances in the polarization ratio of scat- 
tered light. Lettieri et al. [33] measured the diameter of 
droplets in the size range 6 to 12 jam to a resolution of 
0.003 jam using this technique. 

Work is currently in progress to develop a nominal 
0.3 jam particle size standard. In this case, quasielastic 
scattering is beiig used in conjunction with transmission 
electron microscopy with the magnification calibration 
based on the 0.8)5 jam spheres (SRM 1690). 
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APPENDIX 

Preparation and Stability of Polystyrene Sphere Suspensions 



Preparation — The procedure used by Dow Chemical 
Co. to produce monosize polystyrene spheres is known 
as emulsion polymerization. A brief description of the 
general process is provided here, though the specific 
emulsifying agent and chemical initiator used by Dow 
are confidential. A more detailed description of emul- 
sion polymerization is given by Blackley [34]. 

Emulsion polymerization involves four components: 
water, styrene monomer, an emulsifying agent, and a 
free radical initiator. A typical emulsifying agent is so- 
dium lauryl sulfate. The emulsifier will permit the dis- 
persion of the stirred monomer phase into a stable sys- 
tem of micrometer-size droplets. At a sufficiently high 
concentration of emulsifier, hydrated aggregates called 
micelles, made up of 50 to 100 molecules of emulsifier, 
will form. The concentration of micelles is approxi- 
mately six orders of magnitude greater than that of the 
monomer droplets. 

The polymerization reaction begins when the free 
radical initiator is heated to sufficiently high tem- 
perature to dissociate into radical ions. This is around 50 
°C for persulfate anions. The sulfate radical anions at- 
tack the styrene dissolved in the aqueous phase to form 
new ion radicals. The polymer chain begins growing by 
the addition of more styrene molecules to the radicals in 
the aqueous phase. The polymer-containing radicals 
soon enter the micelles where propagation of the chain 
is continued by attack on the monomer within the mi- 
celle. The growth of the polymer chain is supported by 
the rapid diffusion of the styrene from the droplets into 
the growing micelles. This growth is terminated by re- 
action with a second free radical. At this point the chain 
polymer may have reached a molecular weight on the 
order of 10 6 . Later the chain growth is reinitiated by 



another radical, and this process will continue until most 
of the styrene is consumed. 

While the polymer growth is a chain process, the 
intermolecular forces and surface tension lead to a 
spherical shape analogous to a ball of string. The nar- 
rowness of the size distribution results from the uni- 
formity of styrene concentration throughout the reactor 
and the growth time being long compared to the ini- 
tiation time. The concentration of the initiator and emul- 
sifier are also important for obtaining a narrow size dis- 
tribution. Too much emulsifier can lead to new particle 
formation. This is especially true when small poly- 
styrene spheres are used as seed particles for growing 
larger ones. The growth of 1-fim-diameter spheres for 
SRM 1690 required the use of seed particles, and there 
was evidence of a new generation of small particles. 



Stability — The electrostatic charge on the surface of 
the particles, created by the outward pointing ionic end- 
groups, is crucial to the stability of the suspension. 
Changes in the surface layer can lead to agglomeration 
of the spheres. Dilution and the addition of electrolyte, 
required for some of the measurement techniques, can 
affect the stability of the suspension. The Verwey- 
Overbeek theory [35] describes the stability of colloidal 
suspensions in terms of the electrostatic forces of re- 
pulsion between colloidal particles and London-van der 
Waals forces of attraction. The electrostatic repulsion 
arises from the surface charge; for example, from the 
adsorbed emulsifier ions. The repulsion is greatest at the 
point of contact between the electrical double layers on 
the two particles. This is on the order of a nanometer 
from the polymeric surface, since the ionic portion of 
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the emulsifier is directed away from the sphere. At 
greater distances between the particles' centers, the re- 
pulsion drops off exponentially. The magnitude of the 
repulsion depends on the surface particle charge and 
electrolyte concentration. The London-van der Waals 
forces are greatest close to the particle surface (a few 
tenths of a nanometer from the polymeric surface) and 
are not strongly affected by the surface charge of elec- 
trolyte concentration. Of course, there are very strong 
repulsive forces as the polymeric molecules of two 
spheres come in contact. 

The sum of these three forces leads to a potential with 
a minimum at short distances, a repulsive barrier due to 
electrostatic repulsion at greater distances, and possibly 
a second minimum at relatively large distances, resulting 
from the longer range of the London-van der Waals 
forces. If the particle has sufficient energy to penetrate 
the electrostatic barrier of another particle, the two par- 
ticles will stay together in an aggregate. An increase in 
the electrolyte concentration or a decrease in the sur- 
face charge density caused by diluting the suspension 
with electrolyte will lower the repulsive barrier and 
may lead to agglomeration. Such an effect was a con- 
cern in the electrical-sensing zone counter measurments 
as discussed in subsection 5.2, since the particle sus- 
pension was diluted with electrolyte. 
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